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The star formation history of RCW 36 



L. E. Ellerbroek 1 , L. Kaper 1 , A. Bik 2 , K. M. Maaskant 1 , L. Podio 3 ' 4 

Abstract. Recent studies of massive-star forming regions indicate that they can contain mul- 
tiple generations of young stars. These observations suggest that star formation in these regions 
is sequential and/or triggered by a previous generation of (massive) stars. Here we present new 
observations of the star forming region RCW 36 in the Vela Molecular Ridge, hosting a young 
cluster of massive stars embedded in a molecular cloud complex. In the periphery of the clus- 
ter several young stellar objects (YSOs) are detected which produce bipolar jets (HH 1042 and 
HH 1043) demonstrating that these objects are still actively accreting. The VLT/X-shooter spec- 
trum of the jet structure of HH 1042 provides detailed information on the physical conditions 
and kinematical properties of the jet plasma. From this information the YSO's accretion history 
can be derived. Combining the photometric and spectroscopic observations of RCW 36 gives 
insight into the formation process of individual stars and the star formation history of this young 
massive-star forming region. 



1. Introduction 

Important questions in the hitherto poorly understood process of (massive) star formation are: 
(i) How are the most massive stars formed? (ii) Which processes trigger star formation? We 
approach these questions by studying young stars and their environment during or directly after 
the time of formation. Near-infrared imaging surveys provide a census of the embedded young 
stellar populations. Follow-up near-infrared spectroscopy provides a window on the photo- 
spheres of young massive stars, as the extinction is orders of magnitude less than at UV and 
optical wavelengths, and thermal dust emission has not set in. This then results in accurate 
spectral classification and allows for a detailed investigation of the signatures of star formation 
(e.g. disks and jets). 

Massive stars (M > 8 M Q ) are rare, their formation timescale is short and their birth 
places are obscured by dust, leading to high extinction (A v ~ 10 - 100 mag). These factors 



make it difficult to observe massive stars in formation; see e.g. Zinnecker & Yorke (2007) 
for a review. Massive stars are expected to form in a similar way to low-mass stars, with 
the difference that the strong UV radiation field of the contracting massive star will hinder the 
accretion process. The first detectable phase is thought to be the ultra-compact H n region, when 
the young massive star starts to ionize the surrounding medium. Whether or not the central 
star(s) embedded in an UCHu are still contracting to the main sequence and carry the signatures 



of the formation process remains to be demonstrated. Theoretical models by |Hosokawa et al. 



( 2010 ) predict that massive stars can form through disk accretion, though require high accretion 
rates (M > 10" 4 M yr" 1 ). During this rapid accretion process, the radius of the star will 
increase dramatically due to the high entropy deposited on the st ellar surface. At later sta ges, the 
protostar begins to contract towards the main sequence. Thus, |Hosokawa et al.| ( |2010) predict 
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Figure 1. Left: NTT/SOFI narrowband image of RCW 36 (red: Bry, gree n: H 2 , blue: Pa/?). 
Overplotted contours are BLAST 500 iim flux, tracing the molecular cloud (Netterfield et al. 
2009 ). X- shooter targets are labeled. The region within the dashed square is displayed in 
Fig. [4] Right: Hertzsprung-Russell diagram for RCW 36, with an adopted distance of 0.7 kpc. 
Classified stars which were also observed with X-shooter are plotted with red symbols. Dotted 
lines are pre-main sequence (PMS) tracks (5, 4, 3, 2.5, 2, 1.5, 1 and 0.5 M , from top to 
bottom) and the dash-dotted line is the 1 Myr PMS isochrone ( Siess et al. 2000 ). The intrinsic 
near-infrared excess of 408 was subtracted from its K-band flux. 



the existence of bloated high-mass protostellar objects, which lack detectable H n regions. An 
example of such an o bject is B275 in M17, a n 8 M© pre-main sequence (PMS) star with the 
dimensions of a giant ( [Qchsendorf et al.|2011| ). 

Most massive stars form in clusters (e.g., Lada & Lada 2003 Zi nnecker & Yorke|2007 ). 
Both during and after their formation, they will have impact on their surroundings. In some 
well-studied massive star-forming regions, the massive stars are centrally located and appear 



to have formed later than an older generation of intermediate-mass stars (e.g., Feigelson & 



Towns ley 2008; Bik et al.|20 10; Maaskant et al.|2011) . Sequential star formation is not a new 
phenomenon (see e.g. |Blaauw|1964||1991| ). An unanswered question is whether massive stars 
play an important role in this sequence; whether a new generation is triggered by their emergent 
radiation field, stellar winds and supernovae, or perhaps that massive stars are the last to form 
in a cluster. Several scenarios exi st for triggered sta r formation; massive stars are expected to 
play a crucial role in this process ( [Andre et al.|20Q9j ). 

In order to understand this process in more detail, a complete census of the outcome of 
star formation is required. Regions located relatively nearby are best suited for this purpose. 
However, nearby star- forming regions contain few or no massive stars (Taurus- Auriga) or their 
star formation process is already in an advanced stage (Or ion), whil e most massive-star formi ng 
regions are distant (> 2 kpc, e.g. M17, |Hanson et al.|1997] 30 Dor, |Walborn & Blades|1997| ). 

In this paper, we present a near-infrared imaging and spectroscopy survey of the massive 
star-forming region RCW 36. The optical to near-infared spectrograph X-shooter, mounted on 
the ESO Very Large Telescope, is used to probe the youngest and the most massive objects in 
this cluster in order to obtain information about their photospheres and circumstellar material. 
The region contains two jet sources; their spectra contain a fossil record of their recent outflow 
history. In section [2] the star forming region, along with its infrared-brightest stars, is intro- 
duced. Section [3] then elaborates on the two jet sources (HH 1042 and HH 1043). The physical 
conditions of HH 1042, along with a simulation of its kinematics with a simple ballistic model, 
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Figure 2. Spectral energy distributions of the X-shooter sample of RCW 36. The spectra are 
normalized to their mean K-band flux. Overplotted in grey are the dereddened spectra of the 
four classified stars. Dashed lines are Kurucz models corresponding to the spectral type, scaled 
on the dereddened flux. Note the intrinsic infrared excess redward of 1000 nm in 408. 



are presented in section |3.2| The results are put in the context of sequential star formation in 
section |4] 
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Figure 3. Normalized spectra in the Ca n triplet region of the stars in RCW 36 observed with 
VLT/X- shooter. Note the absence of the diffuse interstellar band in the younger objects. 



2. The star-forming region RCW 36 



A good opportunity to study the sequence of star formation is presen ted by the cluster R CW 36, 
part of the Vela C molecular complex, at a distance of 0.7 ± 0.2 kpc (Lisea uet al.|1992| ). Photo- 
metric data from NTT/SOFI and Sp itzer/IRAC, integral field spectroscopy from VLT/SINFONI, 
spectrosco py from VLT/X-shooter ([Ellerbroek et al.|2Q12aj ) and Herschel/PACS and SPIRE ob- 
servations (Hill et al. 2011; Giannini et al. 2012 ) reveal a rich stellar population and cloud 
environment (Fig.jJJ left). The center of the star forming region hosts a young cluster including 
a number of OB stars (|Bik et aL||2005|), while also numerous Young Stellar Objects (YSOs) 



are identified, two of which have strong CO (v = 2 - 0) emission produced by a dense disk 
( |Bik & Thil|2004| |Wheelwright et al.||201 |). Both YSOs produ ce extensive jets with atomic 
gas velocities of several hundred km s 1 ( Ellerbroek et al.|201l| , indicating that they are still 
actively accreting. The region suffers from a high optical exinction due to the local prese nce of 
interstellar dust. By dereddening all sources to the JHK colors of classical T Tauri stars ( Meyer 
et al.|1997]) we find an aver age value of Ay = 10.5 mag, somewhat higher than A v = 8.1 mag 
found by Baba et aL] ( |2004| ) by the same method. Possible explanations for this difference are 
the better observing conditions and higher completeness (up to J =19) of the dataset used in the 
current study. 
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The spectral type and luminosity class of 58 sources were determined by their H- and 
^-band spectra. Four early-type stars were observed with VLT/X-shooter (300 - 2500 nm) 
and classified by their optical spectra. The results of this classification are contained in the 
Hertzsprung-Russell diagram (HRD, Fig.[T] right). A pre-main sequence population is scattered 
around the 1 Myr isochrone, while the two central OB stars (413, 09.5V and 462, 09V) are 
already on the main sequence. The distribution of the stellar masses according to the evolution- 
ary tr acks is consistent with the Salpeter initial mass function for stars in the Galaxy ( Salpeter 
[1955]). 

In Figs. [2] and [3] spectral energy distributions (SED) and normalized spectra are displayed 
of the most significant sources in the cluster. The two O stars 413 and 462 are almost featureless, 
save for diffuse interstellar bands (DIBs). These DIBs are also seen in the other two classified 
stars, while they are absent in the YSO spectra. 

The two most massive PMS stars are 179 (B8) and 408 (AO), on evolutionary tracks of 4 
and 5 M Q , respectively. The Spitzer observations of 179 show no infrared excess at all, nor does 
it exhibit any emission lines. We conclude it is a bloated PMS star, not yet on the main sequence. 
Object 408 has an infrare d excess red ward of 1000 nm (see Fig. [2]) with colors characteristic 
of a Lada class II source ( |Lada|19 87). It also exhibits H i emission lines superposed on broad 
absorption lines. Towards lower transitions the lines become optically thick and have P-Cygni 
profiles indicative of a wind origin. The Ca n triplet profiles have an extended red wing, also a 
wind diagnostic. 

The three other sources in Fig. [2] and [3] are identified as YSOs by their infrared excess and 
pure emission line spectra indicative of disks and outflows. Object 506 is an embedded class 
I source; its H i emission lines have extended red wings. Object 292 is an actively accreting 
system, its photosphere veiled by continuum emission from the accretion disk and its spectrum 
dominated by emission lines showing signatures of Keplerian rotation and outflow. The spec- 
trum exhibits many H i, He i and metallic lines, and the 2290 nm CO bandheads are in emission. 
Also, its many forbidden lines are associated with a (disk) wind or jet, which is also visible in 
the VLT/SINFO NI [Fe n] image (HH 104 2, Fig. |4j. For an extensive discussion on the nature 



of object 292, see|EUerbroek et al. (2011 



Object 480 has a spectrum (not shown in Fig. [3]) dominated by emission lines from both 
the nebula and a jet. Only the Ca n triplet, CO bandheads and a weak continuum redward of the 
H band can be related to the central object. Like 292, it powers a jet (HH 1043) visible in many 
H i and forbidden emission lines. 

From top to bottom, the spectra in Fig. [3] show a tentative sequence of star formation, 
from embedded objects with disks and outflows, via PMS stars with and without disks and 
winds, to main sequence stars with photospheric spectra void of any emission lines. Also note 
the increasing strength of the DIB feature at 862 nm along this sequence. The fact that these 
objects are all found in the same cluster rises the question whether this is evidence for sequential 
or even triggered star formation. No direct evidence for binarity is found in these data. In the 
remainder of this paper, we focus on two individual objects, namely the jets HH 1042 and HH 
1043, powered by 292 and 480, respectively. 



3. HH 1042 and HH 1043 

Jets are a universal signature of angular momentum transport in accreting systems. The ob- 
servational signature of jets emitted by YSO are the so-called Herbig Haro objects. These are 
ionized nebulae produced by the shocks occurring along the jet and become optically visible in 
the late stages of the star formation process, when the environment has been cleared by the jet 
itself. Stellar jets are thought to be powered by a stellar or disk wind, carrying away material 
from the disk along magnetic fieldlines by the centrifugal force. The confinement of the plasma 
along magnetic field lines results in a collimated flow (Blandford & Payne 1982 ). 

The gas heated behind the shock fronts produce emission lines which can be used to re- 
trieve the gas physical conditions and kinematics. For example, the analysis of the velocity 
variation along the flow allows us to investigate the formation of shocks. These can have mul- 
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Figure 4. Left: Composite linemap of the dashed region indicated in Fig. [I] Colors trace the 
[Fe n] il644 nm and H 2 /12121 nm lines as observed with VLT/SINFONI, and the Spitzer/IRAC 
8 fim continuum. The [Fe n] emission traces the jets. Right: Composite linemap of the X- 
shooter spectrum of HH 1042. Colors trace three emission lines ([O i], [N n] and [Fe n]); the 
fluxscale is normalized to the lcr noise level. The central source continuum at 0" is subtracted 
using a Gaussian fit. Note the difference between the [O i] and [Fe n] emission. The emission 
at zero velocity is of nebular origin. 



tiple origins, e.g. a variable outflow rate at the base or interaction with the ambient nebula . In 
jj |3.1| we briefly present the characteristics of the two Herbig-Haro objects in RCW 36. In § 3.2 



we reconstruct the likely accretion history from the kinematics of the gas in one of these jets. 

|Ellerbroek et al.| ( [2011| ) report the detection of two jets in RCW 36 in the [Fe n] 1644 nm 
and Bry images obtained with VLT/SINFONI. They have been classified as a Herbig-Haro 
objects HH 1042 and HH 1043 (B. Reipurth, private communication). The left panel of Fig. [4] 
shows the location of the jets in the star forming region. Note that they are located in the 
periphery of the cluster, along the filaments traced by Bry and H2 emission (Fig. [I] left). A lack 
of mid-infrared continuum emission is detected around HH 1042, which suggests that this jet 
is protruding from the molecular cloud. The northwestern (redshifted, or 'red') lobe fades into 
an ISM structure, while the southeastern ('blue') lobe seems to be colliding with the ISM. It is 
not clear whether the patch of [Fe 11] emission 8" southeast of this location is related to the jet. 
HH 1043 is more embedded; both mid-infrared and H2 emission is detected near with the jet. 
The western lobe terminates in a bow shock feature dominated by [Fe 11] emission. 

Optical to near-infrared, medium-resolution (25 < Av < 50 km s" 1 ) spectra of the two jets 
were obtained with VLT/X- shooter, with the slit aligned along the jet axis.. In the next sections 
the X-shooter spectrum of HH 1042 is presented and discussed. 

3.1. Physical and kinematical structure of the HH 1042 jet 

The VLT/X-shooter spectrum of HH 1042 shows emission in more than 90 emission lines. 
Line velocity variations on an arcsecond scale suggests the interaction of gas components with 
different velocities. Fig. [4] (right panel) shows a composite position velocity (PV) diagram of 
three emission lines which trace different gas components in the flow. Note that the emission 
lines in the red lobe are suffering from extinction due to the molecular cloud. 

The [N 11] line, tracing a medium with a high degree of ionization, dominates at km s" 1 . 
This component is likely background emission from the H 11 region and is unrelated to the jet. Fe 
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is strongly depleted in the ISM and released in gas-phase at the jet launch and/or in the shocks 
occurring along the jet. Hence, the [Fe n] line is exclusively detected in the jet and reveals its 
kinematic structure. We observe a varying velocity along the flow; at -2" and +6" shocks form 
where fast material is 'ploughing' into slower material. These regions are traced by [O i] and 
[N n] emission. At the base of the jet, before strong ionization is produced by the shocks, we 
detect bright [O i] emission, possibly tracing the disk wind origin of the jet. Its velocity at the 
base coincides with the blueshifted absorption component of the Ca n triplet lines (see Fig. [3]), 
strongly suggesting a connection between the jet and a disk wind, although these spectral lines 
may trace different parts of the flow. 

3.2. The accretion history of HH 1042 

We now present the analysis of the kinematics (positions and velocities) of the brightest jet 



tracer in the spectrum ([Fe n] A 1644 nm). We adopt a method first used by Rag a et aL] ( |1990| ) to 
determine the outflow history of the jet. A concise model is conceived in order to simulate the 
flow over the past 900 years. 

We start by measuring the positions on the sky x t and radial velocities v r of the local emis- 
sion maxima, or 'knots' (see Fig. [5] lower right panel). Assuming that the jet is perpendicular 
to the disk, and ado pting the value / = 18° (the inclination of the jet axis with respect to the line 
of sight) derived by Wheelwright et al.| ( 2010| ), we can then calculate for each knot a timescale, 



which can be interpreted as the time at which the knot was launched, given no collisions have 
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Figure 5. Model simulation assuming a flow with parameters specified in the text. The 
observed [Fe n] A\ 644 emission position-velocity diagram (bottom-right panel) is compared 
with the result from the model simulation. The adopted outflow rate and its parameters are 
specified in the text. The dashed lines indicate the knot positions in the observed and simulated 
position- velocity diagrams. 
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occured in the flow: ^ 

^launch — " : (1) 

v r tan i 

We then fit a periodic outflow rate at the jet base to the resulting values of v = v r / cos i and 

^launch- 



v(0 = 170 km s _i + 42 km s _i sin 2tt (2) 

\ 280 yr / 

We furthermore note that the knots in the red lobe are fit better by this function if it is offset 
by a phase difference of n/2. Also, an additional mode with a period of 32 yr is introduced 
to produce the substructure along the flow. We use this function as an input for a ballistic 
model, where no collisions occur along the flow, assume a constant mass loss rate, and run a 
simulation for t = [-1000, 0] yr (Fig. [5]). The global emission profile of the [Fe n] line can be 
reproduced reasonably well. It should be noted that the colors in the simulation images in Fig. [5] 
trace the material density, while the observed quantity is line emission. The derived velocities 
and timescales, as well as the apparent phase shift between the two lobes, can put constraints 
on theoretical models fo r the jet launching mechanism. For a detailed discussion we refer to 
Elle rbroeketaLlp012b| ). 



4. Discussion 

The stellar population of RCW 36 consists of a large number of low-mass PMS stars (< 2M©), 
a couple of intermediate-mass PMS stars (2-8 M ) and two high-mass stars (~ 20 M Q ) already 
on the main sequence. In addition to this, at least three YSOs are present, of which the spectral 
type cannot be determined; based on their spectral char acteristics (CO, H i, He i) their central 
objects are luminous intermediate- to high-mass stars (Ellerb roek et al.|2()TT| ). Two of these are 
associated with HH outflows, indicating that accretion is still ongoing. 

The simultaneous presence of low-mass PMS stars and intermediate-mass YSOs suggests 
that these belong to different generations of star formation. The low- to intermediate-mass PMS 
stars seem to have formed about 1 Myr ago. The O stars are close to the ZAMS, implying they 
are of an equal age or youn ger (< 1 Myr), if thei r PMS phase is indeed as rapid as suggested 
by observations (~ 10 5 yr, |Mottram et al.||201l] ). Shortly after having formed they possibly 
triggered the formation of the YSOs in the periphery of the cluster. This is also suggested by 
the presence of a shock wave (the filament traced by Bry and H2 emission, Fig. [T] left). The 
dynamical timescale associated with the distance between the O stars and the shock wave (0.2 
pc), assuming it has propagated with a typical soundspeed of a few km s" 1 , is of the order of 
10 4 yr. This, and the young age of the YSOs, is consistent with the YSOs having formed after 
the massive stars. However, this scenario cannot be confirmed, mainly because the age of the 
O stars is poorly constrained due to the uncertainty in the distance to the cluster (i.e. the precise 

location of the ZAMS). 

Recent studies ( [Feigelson & Townsley|2008l |Bik et al.|2010|[Wahg et al.|20li]|Maaskant 



et al. 



2011) have found similar age spreads in massive-star forming regions, with the most 



massive stars being the youngest in the cluster. Whether these observations really suggest an 
evolutionary sequence, and whether or not high-mass stars form after lo w-mass stars, remains 
a matter of debate; for a discussion see e.g. Zinnecker & Beuther (2008 ). 



VLT/X- shooter is a powerful instrument to study star formation, since both the stellar 
photosphere and the circumstellar material are probed with a single observation. This is also 
illustrated by the results of the kinematic study of HH 1042. The observed velocity structure in 
the jet can be reproduced reasonably well by assuming a ballistic flow with a periodic outflow 
rate. Incorporating relevant line excitation mechanisms in the model might improve the match 
between the observed intensities of the knots and the results from the simulation. 

The observed jet velocity structure is consistent with the escape velocity from a launch 
region within a few stellar radii of an intermediate-mass star. The model parameters that best 



reproduce the observed emission in the position- velocity diagram are comparable with those 
seen in other jets (e.g. Velazque z et al. [2004{ Raga et al. 2012 ). The variations in velocity 
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might then be caused by episodes of accretion, being of the order of the viscous timescale in 
the accretion disk. Another explanation would be that the observed velocity oscillation is just a 
projection of the global movement of the system: a binary orbit, or jet precession. 

The phase shift in velocity between the red and blue lobe would follow naturally from a 
precession, although the phase shift should be n rather than 7r/2, which is within the uncertainty 
limit of the velocity of the central source. Also, the precession angle cannot be larger than -10° 
due to the observed collimation of the jet. Another explanation for an asymmetry between 
the two lobes would be the different environmental conditions on either side of the system 
flMelnikov et al.|2009t|Podio et al.|2011| ). 

Questions. R. Oudmaijer: On what timescale can we detect proper motions of the HH 1042 ? 

With the limited means of performing astrometry on the SINFONI and X-shooter observa- 
tions, which were taken 4 years apart, we can at best detect proper motions in the order of 1". 
These are not found. This corroborates the conclusion of the kinematics analysis, yet it does 
not put a firm constraint on the inclination. Based on our results we predict proper motions in 
the order of 0.1" on a timescale of ~ 5 yr. 

C. Martayan: Do you have an estimate of the inclination and spatial extent ofHH 1042? 

Adopting a distance of 0.7 kpc, the projected length of the jet is 1.2 x 10 4 AU. Lacking 
proper motion measurem ents, we have derived the jet inclination from the disk inclination (18°, 
Wheelwright et al. 2010) by assuming the jet is collimated and perpendicular to the disk equato- 
rial plane. Using this value the deprojected length is 3.9 x 10 4 AU. However, the adapted value 
for the incli nation might be an underes timate due to a possible overestimated mass of the central 
star used by Wheelwright et al. ( 2010 ). Note that the quoted value is the length of the optically 
visible part of the jet; its true extent may be much larger and visible at longer wavelengths. 

References 

Andre, P., Basu, S., & Inutsuka, S. 2009, The formation and evolution of prestellar cores (Cambridge 
University Press), 254 

Baba, D., Nagata, T, Nagayama, T., Nagashima, C, Kato, D., Kurita, M., Sato, S., Nakajima, Y., Tamura, 

M., Nakaya, H., & Sugitani, K. 2004, ApJ, 614, 818 
Bik, A., Kaper, L., Hanson, M. M., & Smits, M. 2005, A&A, 440, 121 

Bik, A., Puga, E., Waters, L. B. F. M., Horrobin, M., Henning, T., Vasyunina, T, Beuther, H., Linz, H., 
Kaper, L., van den Ancker, M., Lenorzer, A., Churchwell, E., Kurtz, S., Kouwenhoven, M. B. N., 
Stolte, A., de Koter, A., Thi, W. E, Comeron, E, & Waelkens, C. 2010, ApJ, 713, 883 

Bik, A., & Thi, W. E 2004, A&A, 427, L13 

Blaauw, A. 1964, in The Galaxy and the Magellanic Clouds, edited by E J. Kerr, vol. 20 of IAU Sympo- 
sium, 50 

— 1991, in NATO ASIC Proc. 342: The Physics of Star Formation and Early Stellar Evolution, edited by 

C. J. Lada & N. D. Kylafis, 125 
Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883 

Ellerbroek, L. E., Bik, A., Maaskant, K. M., Kaper, L., Paalvast, M., & Balog, Z. 2012a, in prep. 
Ellerbroek, L. E., Kaper, L., Bik, A., de Koter, A., Horrobin, M., Puga, E., Sana, H., & Waters, L. B. E M. 
2011, ApJ, 732, L9 

Ellerbroek, L. E., Podio, L., Kaper, L., Sana, H., & de Koter, A. 2012b, in prep. 
Feigelson, E. D., & Townsley, L. K. 2008, ApJ, 673, 354 

Giannini, T, Elia, D., Lorenzetti, D., Molinari, S., Motte, E, Schisano, E., Pezzuto, S., Pestalozzi, M., 
di Giorgio, A. M., Andre, P., Hill, T, Benedettini, M., Bontemps, S., di Francesco, J., Fallscheer, 
C, Hennemann, M., Kirk, J., Minier, V., Nguyn Lu'O'Ng, Q., Polychroni, D., Rygl, K. L. J., 
Saraceno, P., Schneider, N., Spinoglio, L., Testi, L., Ward-Thompson, D., & White, G. J. 2012, 
A&A,539,A156 

Hanson, M. M., Howarth, I. D., & Conti, P. S. 1997, ApJ, 489, 698 

Hill, T, Motte, E, Didelon, P., Bontemps, S., Minier, V., Hennemann, M., Schneider, N., Andre, P., 
Men'shchikov, A., Anderson, L. D., Arzoumanian, D., Bernard, J.-P, di Francesco, J., Elia, D., 
Giannini, T, Griffin, M. J., Konyves, V., Kirk, J., Marston, A. P., Martin, P. G., Molinari, S., 
Nguyen Luong, Q., Peretto, N., Pezzuto, S., Roussel, H., Sauvage, M., Sousbie, T, Testi, L., 
Ward-Thompson, D., White, G. J., Wilson, C. D., & Zavagno, A. 2011, A&A, 533, A94 



10 



Hosokawa, T., Yorke, H. W, & Omukai, K. 2010, ApJ, 721, 478 

Lada, C. J. 1987, in Star Forming Regions, edited by M. Peimbert & J. Jugaku, vol. 115 of IAU Sympo- 
sium, 1 

Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57 

Liseau, R., Lorenzetti, D., Nisini, B., Spinoglio, L., & Moneti, A. 1992, A&A, 265, 577 

Maaskant, K. M., Bik, A., Waters, L. B. F. M., Kaper, L., Henning, T., Puga, E., Horrobin, M., & Kainu- 

lainen, J. 2011, A&A, 531, A27 
Melnikov, S. Y., Eisloffel, J., Bacciotti, F, Woitas, J., & Ray, T. P. 2009, A&A, 506, 763 
Meyer, M. R., Calvet, N., & Hillenbrand, L. A. 1997, AJ, 114, 288 

Mottram, J. C, Hoare, M. G., Davies, B., Lumsden, S. L., Oudmaijer, R. D., Urquhart, J. S., Moore, 
T. J. T., Cooper, H. D. B., & Stead, J. J. 201 1, ApJ, 730, L33 

Netterfield, C. B., Ade, P. A. R., Bock, J. J., Chapin, E. L., Devlin, M. J., Griffin, M., Gundersen, J. O., 
Halpern, M., Hargrave, P. C, Hughes, D. H., Klein, J., Marsden, G., Martin, P. G., Mauskopf, R, 
Olmi, L., Pascale, E., Patanchon, G., Rex, M., Roy, A., Scott, D., Semisch, C, Thomas, N., Truch, 
M. D. P., Tucker, C, Tucker, G. S., Viero, M. P, & Wiebe, D. V. 2009, ApJ, 707, 1824 

Ochsendorf, B. B., Ellerbroek, L. E., Chini, R., Hartoog, O. E., Hoffmeister, V., Waters, L. B. F M., & 
Kaper, L. 201 1, A&A, 536, LI 

Podio, L., Eisloffel, J., Melnikov, S., Hodapp, K. W, & Bacciotti, F. 2011, A&A, 527, A13 

Raga, A. C, Binette, L., Canto, J., & Calvet, N. 1990, ApJ, 364, 601 

Raga, A. C, Rodriguez-Gonzalez, A., Noriega-Crespo, A., & Esquivel, A. 2012, ApJ, 744, L12 
Salpeter, E. E. 1955, ApJ, 121, 161 

Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 
Velazquez, P F, Riera, A., & Raga, A. C. 2004, A&A, 419, 991 
Walborn, N. R., & Blades, J. C. 1997, ApJS, 112, 457 

Wang, Y, Beuther, H., Bik, A., Vasyunina, T., Jiang, Z., Puga, E., Linz, H., Rodon, J. A., Henning, T., & 

Tamura, M. 2011, A&A, 527, A32 
Wheelwright, H. E., Oudmaijer, R. D., de Wit, W J., Hoare, M. G., Lumsden, S. L., & Urquhart, J. S. 

2010, MNRAS, 408, 1840 

Zinnecker, H., & Beuther, H. 2008, in Massive Star Formation: Observations Confront Theory, edited by 
H. Beuther, H. Linz, & T. Henning, vol. 387 of Astronomical Society of the Pacific Conference 
Series, 438 

Zinnecker, H., & Yorke, H. W. 2007, ARA&A, 45, 481 



